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ABSTRACT 

Typical gamma-ray burst spectra are characterized by a spectral break, E p , which for bright BATSE 
bursts is found to be narrowly clustered around 300 keV. Recently identified X-ray flashes, which may 
account for a significant portion of the whole GRB population, seem to extend the E p distribution to 
a broader range below 40 keV. Positive correlations among E p and some other observed parameters 
have been noticed. On the other hand, within the cosmological fireball model, the issues concerning 
the dominant energy ingredient of the fireball as well as the location of the GRB emission site are 
still unsettled, leading to several variants of the fireball model. Here we analyze these models within a 
unified framework, and critically reexamine the E p predictions in the various model variants. Attention 
is focused on the predictions of the narrowness of the E p distribution in different models, and the 
correlations among E p and some other measurable observables. These model properties may be tested 
against the current and upcoming GRB data, through which the nature of the fireball as well as the 
mechanism and site of the GRB emission will be identified. In view of the current data, various models 
are appraised through a simple Monte-Carlo simulation, and a tentative discussion about the possible 
nature of X-ray flashes is presented. 

Subject headings: gamma rays: bursts - radiation mechanisms: non-thermal - shock waves - stars: 
magnetic fields 



1. INTRODUCTION 

Gamma-ray burst (GRB) sub-MeV spectra contain 
most of the prompt information available from these mys- 
terious sources. Progress in understanding the origin of 
such spectra, however, has been slow. A typical GRB 
spectrum is non-thermal, and can be fitted by a so-called 
Band-function (Band et al. 1993) with three parameters: 
a low energy power law photon index a, a high energy 
power law photon index /3, and the spectral break energy 
E p which defines the smooth transition between the two 
power laws. Since generally (3 < —2, E p is related to 
the peak of the vF v spectrum, and therefore is also called 
"E-peak" . Within the framework of the commonly con- 
sidered cosmological fireball model, there are a number 
of variants (invoking different fireball contents, different 
emission sites, or different emission mechanisms) proposed 
to explain the GRB data. An important quantity which 
characterizes the GRB is E p , which for 156 bright BATSE 
bursts in the 4B sample is found to be narrowly clustered 
around 300 keV, with a log-normal distribution with full 
width at half maximum of less than a decade (Preece et al. 
1998; 2000). Theoretically, the value of E p is expected to 
be correlated with some other observational parameters, 
but such correlations could vary significantly in different 
models. These provide important criteria to evaluate the 
correctness of the models by comparing such correlations 
found in the data. 

This topic has received heightened attention in view of 
two recent developments. On the one hand, a new category 
of X-ray transients, known as X-ray flashes (XRFs) has 
been identified (Heise et al. 2001; 2002, in preparation). 
The XRFs resemble normal GRBs in many respects, with 
the novelty that the peak energies are distributed from 
100 keV to below 40 keV. These objects appear to form a 



natural extension of the GRB population in the softer and 
fainter regime (Kippen et al. 2001; 2002), which widens 
the narrow E p distribution found in the BATSE data, and 
it is estimated that XRFs represent a large portion (e.g. 
^1/3) of the whole GRB population. It is therefore in- 
teresting to see how present theories can accommodate 
this new category of GRBs. On the other hand, GRB 
light-curve variabilities (Fenimore & Ramirez-Ruiz 2000; 
Reichart et al. 2001) and spectral lags (Norris et al. 2000; 
Norris 2002) have been proposed as Cepheid-like luminos- 
ity indicators for the long duration (ib>,2 s) GRBs. These 
offer the exciting prospect of a direct relation between ob- 
servable quantities and some of the theoretically most rel- 
evant parameters, such as the wind luminosity L, GRB 
intrinsic durations, emission-site magnetic fields, as well 
as the characteristic of the synchrotron or inverse Comp- 
ton energies (which in many models are linked to the ob- 
served E p ). Empirically, some such correlations have been 
noticed. Lloyd-Ronning & Ramirez-Ruiz (2002) discov- 
ered a positive dependence of E p on the GRB variability 
(or luminosity). When combining this with the luminosity 
- variability correlation (Fenimore & Ramirez- Ruiz 2000; 
Reichart et al. 2001), one infers a positive correlation be- 
tween E p and the burst luminosity L. Indeed, such a corre- 
lation has been seen in the BeppoSAX bursts with known 
redshifts (Amati et al. 2002). Thus, it is timely to criti- 
cally revisit the physical E p predicted in various models, 
as a first step towards the goal of constraining or even 
identifying the nature of the fireball as well as the relevant 
emission site and mechanism for the GRB prompt emis- 
sion. This is the purpose of the present paper. In §2.1, 
we present a synthesis of the current GRB model vari- 
ants within a unified framework, and define the parameter 
regimes in which each model variant applies. In the rest of 
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§2, we revisit these model variants, focusing specifically on 
the i? p predictions. In §3, these predictions, as collected in 
Table g, are used to evaluate the current models through 
a simple Monte-Carlo simulation, and the possible nature 
of the X-ray flashes is tentatively discussed. Conclusions 
are drawn in §4. 

2. GRB MODELS AND E P PREDICTIONS 

2.f . GRB model synthesis 

Cosmological fireball models invoke a brief release of 
energy ~ 1 51 — 1 53 ergs within a short duration of time 
~ 10~ 3 — 10 3 s. The fireball must be clean (low baryon 
load) so that after initial acceleration, the fireball is rela- 
tivistic. The non-thermal spectrum requires the emission 
to be optically-thin, so that the emission site should be 
above the photosphere defined by the baryon content. The 
fireball eventually decelerates by the interstellar medium 
when the afterglow starts. The site of the GRB emission 
is therefore limited to the regime 



pie a decreases with increasing radius due to conversion 
of part of the Poynting flux into the kinetic energy, es- 
pecially when magnetic reconnection is operating, though 
the effect is likely not to be large (Lyubarsky & Kirk 2001; 
Contopoulos & Kazanas 2002). In any case, in the follow- 
ing discussions a refers to the specific value at the relevant 
radius in the problems. 

The deceleration radius is defined by the kinetic energy 
of the fireball, Ek — E/(l + a), if the kinetic energy is 
decoupled from the magnetic energy by the time of decel- 
eration 3 , so that 
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where n ext is the density of the interstellar medium, T is 
the bulk Lorentz factor of the fireball before deceleration, 
and hereafter the convention Q n — Q/IQ™ is adopted. The 
photosphere radius depends on a number of factors and has 
been discussed in Meszaros et al. (2002) (see also §2.4). 
For a very clean fireball in which the opacity is defined by 
the thermal pairs rather than baryonic electrons, 

r ph , ± =r (e /Q' ± ) = 1.9 xlO 9 cm ^{X + oY^H^l _ 3 , 

where 0' ± ~ 0.03 is the normalized (in unit m e c 2 ) comov- 
ing temperature of the fireball when the pairs drop out of 
equilibrium (Shemi & Piran 1990; Paczynski 1986; Good- 
man 1986), and 

9o = (k/m e c 2 )(L h /4nr 2 a)^ = 1.9L^(l+ ff )- 1 /^_ 3 

(4) 

(a is the Stefan-Boltzmann constant, hereafter symbols 
like m e , m p , c, e, h, etc. denote fundamental physical 
constants with conventional meanings) is the initial tem- 
perature of the fireball at the radius rn ~ c£„ )Tn = 3 x 
10 cm i„ imj _3 (t Vtm is the minimum variability timescale 
of the central engine, which defines the width of each mini- 
shell to be ct Vim ), and Lh — L/(l + a) is the luminosity 
of the hot component. In the case where the baryon load 
is somewhat larger and the opacity is dominated by bary- 
onic electrons, the photosphere radius is given instead by 
(5) of Meszaros et al. (2002). In any case, equation 
J) marks the radius above which the pairs are no longer 
in equilibrium. 

Another radius of interest is the critical radius where 
the MHD treatment breaks down, i.e., r MHD . A convenient 
way to define r MHD is to require that the local plasma den- 
sity, n = rib + n± which includes both the baryon density 
and the pair density, to be equal to the Goldreich-Julian 
(1969, hereafter GJ) density, rigj, which is the minimum 

1 Strongly magnetized central engines have been invoked in GRB models motivated by their ability to launch collimatcd jets and to avoid 
heavy baryon loading (Usov 1992; Meszaros & Rees 1997b; Wheeler et al. 2000). Also, the spindown luminosity of the rapidly rotating central 
object, whether a black hole or magnetar, can be tapped through its Poynting flux, which in many cases is more powerful than the initial 
thermal luminosity of the fireball, e.g. Lee, Wijers &; Brown 2000; van Putten 2001). 

2 Strictly speaking, such a definition is equivalent to that in the pulsar wind nebula theory only when a is not too large (e.g. a < 10 4 ). This 
is because, as we know, a pure cold component (e.g. the pulsar wind from Crab) also has a small fraction [~ (10 -5 -10 -4 )] of e nerg y stored as 
the kinetic energy of the pairs flowing from the pulsar magnetosphere. In this sense, any extremely large a (e.g., <x c 3 below, eq.JlCj]) no longer 
has the conventional meaning of the Poynting flux-to-kinetic energy ratio. 

3 In the MHD regime, the kinetic energy and the magnetic energy are coupled, so the deceleration radius is still defined by the total energy 
E of the fireball. Since we are interested in the comparison between r,i £c and r MHD (eq. [p|, at which both energy components are decoupled), 
it is more relevant to define r,j ec with Ex- 



where r p h (e.g. eq.||, seeeq.[5] of Meszaros et al. 2002 for 
a more general treatment) and r^ec (eq.[Q) are the bary- 
onic photosphere radius and the deceleration radius, re- 
spectively. The non-thermal nature of GRB emission also 
requires that the emission energy is not directly coming 
from the hot fireball (which gives a thermal-like spectrum), 
but derives from some other forms, e.g., the kinetic energy 
of the baryon bulk or the magnetic energy of the fireball. 
Generally, the GRB central engine involves a rapidly spin- 
ning and possibly highly magnetized object such as a black 
hole - torus system or a millisecond magnetar 1 . 

The fireball luminosity therefore may be broadly divided 
into two components. The first component, which is the 
component conventionally invoked in the simplest fireballs, 
is initially composed of thermal photons, pairs and a small 
amount of baryons. This component, which we call the 
hot component, essentially stores its energy in the form of 
the kinetic energy of the baryons after initial acceleration, 
apart from (usually) a small fraction of energy leaking at 
the baryonic photosphere. The second component, which 
we call the cold component, is carried by a Poynting-flux 
and low-frequency waves associated with the central en- 
gine spin. For ease of discussion below, we broadly define 
the hot and the cold luminosity components as Lh — Lk 
and L c ~ Lp, respectively, so that the total fireball lumi- 
nosity is L = Lh + L c = Lk + Lp. We also define the ratio 
of the cold-to-hot components as a = L c /Lh — Lp/Lpc 
following the convention of pulsar wind nebula theories 
(e.g. Rees & Gunn 1974; Kennel & Coroniti 1984) 2 . A 
fireball may be then defined as Poynting-flux-dominated 
if <t>1 or as kinetic-energy-dominated if o~<l. In princi- 
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density required for the plasma to be frozen in the mag- 
netic field (Usov 1994) 4 . The GJ density drops as oc r~ 3 
within the light cylinder, r; c = c/f2 (where fi is the angu- 
lar frequency of the central engine), but as oc r~ l beyond 
ric, so that at large distances, the GJ density is ugj = 
{riB s /2irec){R/r lc ) 3 (r lc /r) = {B s R 3 fl 3 /2irec 3 r), where R 
is the radius of the magnetized central engine. Assuming 
that the spin-down energy of the central engine is mainly 
carried by the Poynting flux, one has Lp = £? 2 i? 6 Sl 4 /6c 3 . 
The G J density in the rest frame of the fireball can be then 
expressed in terms of the measurable parameters (here- 
after the primed quantities denote those measured in the 
rest frame of the fireball) 

n' GJ = 1.0 x 10 8 cm" 3 [a/(l + a)] 1 ' 2 L^X^r^ 1 , 

(5) 

where t Vt7n ~ 2ir/Vl ~ 10 3 s, and the value has been 
normalized to r ~ 10 13 cm which is the typical radius 
for the "internal" 7-ray emission (e.g. from the internal 
shocks). The comoving plasma density at the distance r 
is n' = n' b + n' ± . For a continuous wind, the comoving 
baryon density n' b = L k / '47rr 2 'm p c 3 T 2 is 



i' b = 1.8 x 10 13 cm-' 3 L 52 (l + cr)- 1 r^T 2 



(6) 



The pair density is much lower than this. While they are 
in equilibrium, which is below r p h,± (eq.||) and generally 
somewhere below the baryonic photosphere, the pair num- 
ber density is n' ± (th) ~ 4.41 x 10 30 cm- 3 9' 3/2 exp(-e /_1 ) 
(Paczhski 1986). At the radius where <d'± ~ 0.03, this is 
n'±(th) = 7.6 x 10 13 cm.- 3 . Beyond the pair-freeze radius, 
n'± drops as oc r" 2 , so that at the fiducial radius ris — 1) 
one has (cf. Usov 1994) 



n'± ~ 2.8 x 10" cm- 3 L% 2 {1 + <r)- i r%L,- 3 rl3 '■ (?) 

In principle, further annihilation of these pairs is possible. 
In the context of the current problem, the annihilation 
time scale is much longer than the expansion timescale, so 
(ffl) generally applies. While n' ± < n' GJ , the baryon density 
rr b 3> n' GJ at small radii, although it will drop below n' GJ 
at a large enough radius due to the different r-dependence 
of n' b and n' GJ . By requiring n' b ~ n' GJ , one can define the 
radius where the MHD approximation breaks down, i.e., 

r MHD = 1.8 x 10 18 cm Lg 2 [u{\ + a)]- 1 ' 2 ^.^ 1 . (8) 

This radius is usually beyond the deceleration radius r ( i ec 
(so that MHD approximation never breaks) unless 



1nnr 3/4 ,-,-1/2 ,3/2 1/2^-1/2 

a > <t c2 = 190L 52 E 52 t v ' m _ 3 n e l xt T 2 



(9) 



By setting n' b = n' ± , one can also define a critical value 

a = ( 7 C 3=4.1xl0 9 L 5 < m ^3r- 4 , (10) 

above which the electrons associated with baryons are neg- 
ligible. 

Finally, there is another critical value of 



which separates the regimes where strong shocks can or 
cannot develop. Strong shocks are only possible for low 
a flows. For a > o~ c \, the shocks are quite weak, mainly 
because the three speeds on both sides of the shock are 
close to the speed of light (Kennel & Coroniti 1984). 

A fireball may be classified into several sub-categories 
depending on the value of a. Table [l lists the characteris- 
tics of different fireball variants. GRB models can be also 
naturally classified into three categories based upon the 
proposed location of the GRB prompt emission. Within 
each category, several subtypes are defined according to 
the regime of er, as follows. 

1. Internal models, with r p h < r < rdec- In the standard 
scenario the energy input is an unsteady, kinetic-energy- 
dominated wind which undergoes dissipation in internal 
shocks (Rees & Meszaros 1994), which applies for a < a c \. 
For a > cr c i , the source of energy dissipation may be from 
the strong magnetic fields in the fireball. For a > cr C 2, 
the MHD approximation breaks down beyond the radius 
r MHD < r dec, and a global energy dissipation is expected to 
occur. For the extreme case of a > <t C 3, the baryon content 
is negligible, and the process is analogous to a relativistic 
pulsar wind (Usov 1994; Lyutikov & Blackman 2001). For 
a c i < a < o~ C 2, the global MHD approximation applies all 
the way to the deceleration radius. Some local breakdown 
of the MHD condition is required (and possible), probably 
through magnetic reconnection (e.g. Drenkhahn & Spruit 
2002). 

2. External models, with r — rdec- The standard picture 
is that of an external shock model (Rees & Meszaros 1992; 
Meszaros & Rees 1993; Dermer, Chiang & Bottcher 1999). 
A variant of this model invokes the magnetic wind- medium 
interaction in the high-cr regime (Smolsky & Usov 2000). 

3. Innermost models, with r>,r p ^. Although the thermal 
nature of the photosphere emission seems to preclude its 
role as the main GRB emission mechanism, a few bursts 
are known to have quasi-thermal 7-ray spectra. A non- 
thermal character may also appear in the photosphcric 
emission through Comptonization during the emergence of 
the spectrum (Goodman 1986). Also, with a strong mag- 
netic component (not necessarily of very high a), Alfven 
turbulence induced and propagated from the deep layers 
beneath the baryonic photosphere can also lead to a non- 
thermal Compton tail (Thompson 1994; Meszaros & Rees 
2000). Even without Comptonization, both the baryonic 
and the pair photosphere components can provide addi- 
tional components to the optically-thin component, and 
under certain circumstances even play a dominant role, 
which have been suggested to be responsible for the re- 
cently identified X-ray flashes (Meszaros et al. 2002). 

The most natural radiation mechanism in the optically 
thin region (both for the internal and the external models) 
is the synchrotron radiation (or its variants, e.g., jitter ra- 
diation, synchro-Compton radiation, random electric field 
radiation, etc.) and its self-inverse Compton (IC) emis- 
sion, although the acceleration mechanism and the distri- 
bution function of the emitting electrons may vary in dif- 
ferent models. For the baryonic photosphere Comptoniza- 
tion model, the E p is essentially defined by the thermal 
peak although the spectral shape is modified. Below we 
revisit the E p predictions in these model variants. 



a = a cl ~ (0.1 - 1) , (11) 

4 For an alternative but intrinsically similar discussion about the MHD condition, see Spruit, Daigne & Drenkhahn (2001) 
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2.2. Internal models 

Generally, the characteristic synchrotron emission en- 
ergy (which is naturally connected to the E p energy in 
this model) of the electrons with typical energy 7 e m e c 2 
in a relativistic ejecta of bulk Lorentz factor T is E^ yn = 
(4T/3)(3/2)j^(heB' /m e c)sin^(l + z)' 1 , or 5 



1.4 x 10 6 G (Cies,-!^) 1 / 2 ^ 2 ^ 1 ^ 1 , 



(14) 



jpsyn 
P 



2 x 1(T 8 eV (TB')j 2 sin*(l 



(12) 



where B' is the comoving magnetic fields, W is the mean 
pitch angle of the electrons, and z is the source red- 
shift. The factor (TB 1 ) is the Lorentz-boosted mag- 
netic field strength, whose energy density should be pro- 
portional to the fireball energy density in the observer's 
frame as long as the magnetic energy density is a con- 
stant fraction of the total energy density. This gives 
TB' oc U 1 ' 2 oc (L/r 2 ) 1 / 2 oc L x l\- X , where U and L are 
the total energy density and the luminosity of the fire- 
ball. In principle there could be two components that 
contribute to the co-moving magnetic field in the emission 
region. The wind itself usually carries some "primordial" 
magnetic field (B' w ) from the magnetized central engines, 
the lab- frame energy density of which can be estimated as 
Ub = B^/8-k = Lp /A-ncr 2 . The comoving field strength 
is then 

B' w = [2La/(l + a)er 2 ]^ 2 /T 
~ 8.2 x 10 5 G [a/(l + v^LH 2 ^ 1 . (13) 

Although such a large scale wind field is globally organized, 
we assume that in the dissipation regions (internal shocks 
or the reconnection sites) such a field is re-distributed ran- 
domly with comparable strength. Alternatively, magnetic 
fields could be generated in-situ (B' s ) if the energy dissipa- 
tion is through the internal shocks (e.g. Medvedev & Loeb 
1999). Assuming that each proton gains a mean random 
from the internal shocks, the co-moving 



where r ~ T 2 Ao (the internal shock radius) has been 
used. By comparing B' w and B' sl we find that the wind 
component and the shock component of the magnetic 
fields have a similar amplitude for reasonable parame- 
ters, and share the same dependences on the parameters 
such as L52, ri3 and T. This is understandable, since 
both the wind magnetic density and the shock magnetic 
density are assumed to be a constant fraction of the to- 
tal energy density of the fireball (through a and es, re- 
spectively). For typical parameters, the co-moving field 



strength is B' = B' w + B' s = 10' 



G € x iL 5 ' 2 r 



13 ^2 



where 



e xl = 0.82[cr/(l + cr)] 1 / 2 + 1.4(Cies,-i6'p) 1 / 2 for the internal 
shock models, or e x \ = 0.82[cr/(l + a)] 1 / 2 for the magnetic 
field dissipation models. We then have 
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10 8 G e xl Lli 2 r 



52 ' 13 1 



(15) 



/pi iipt 



which removes the apparent T-dependence (although 7*13 
is still dependent on T in the internal shock models). No- 
tice that even if es -C 1 (as shown in the afterglow fits), 
e x i can be still close to unity given a large enough wind 
component (e.g. a > 0.1). Equation (15) is essentially 
model-independent. 

The characteristic Lorentz factor of the emitting elec- 
trons 7 e , however, depends on the poorly understood 
mechanisms of particle acceleration. Conventionally, given 
a typical magnetic field at the acceleration site (eq.fll5|l), 
there are two characteristic 7 e 's from simple scalings. 
One is defined through the simple equipartition argu- 
ment, i.e., to parameterize the random electron energy 
density to be proportional to the comoving thermal en- 
ergy. Due to usually the same scaling laws of both den- 
sities, the typical electron Lorentz factor obtained us- 
ing this method, denoted as ji , generally does not de- 
pend on the radius of emission and particle accelera- 

(2) 

tion. A second estimate of j e , denoted as 7e , is ob- 
tained by limiting the maximum gyro-acceleration at the 
local magnetic field by the synchrotron (and/or IC) cool- 
ing. For the simplest case, the typical acceleration time 
is t' a = 2Trj e (m e c/eB'), and the typical cooling time is 



t' c = (7 e m e c 2 )/[(4/3)7 2 a T c(i?'787r)], 
Thomson cross section. Solving t' a = 
achievable electron Lorentz factor is 



where <jt is the 
t' the maximum 



7< 2 > = (3/2)m eC 2 /[(2 7 r) 1 / 2 e 3 / 2 B' 1 / 2 ] = 4.7 x 10 4 £?' 



energy 

thermal energy density is v! = n'OpiripC 2 , where n! is the 
comoving baryon number density in the shock-heated re- 
gion. The parameter p ~ (7^ — 1) depends on the relative 
Lorentz factor between the shells, 7^ ~ {"fi/jj +7j/7i)/2, 
where 7; and jj are the Lorentz factors of the two col- 
liding shells, respectively. Though uncertain, 9 p is of or- 
der unity for typical parameters, and is not expected to 
be sensitive to the absolute value of T. When a collision 
occur, the slow shell enters the spreading phase due to 
the internal velocity difference within the shell itself. In 
this regime, the density in the shock-heated region evolves 
as n' ~ CM" r/(4/3)7rr 3 m p - C^/(4/3)7rr 3 m p c 2 oc r~ 3 , 
where Mo = Lt v /Tc 2 is the total mass in the colliding 
shell, t v is the typical injection duration of the mini-shells, 
so that Ao ~ ct v is the typical initial shell width, and C 
is the compressive ratio (at least 7 for strong shocks). If 
the shock-generated magnetic fields reach an equipartition 
fraction cb of the total thermal energy in the shock heated 
region, the co-moving field strength in the forward shocked 
region (in the slow shell) is then 

B' s = {8TTe B n'9 p mpC 2 ) 1/2 

5 By adopting such a notation, one has implicitly assumed that there is relativistic emission beamed towards the line of sight. If the rela- 
tivistic ejecta is in a jet with a sharp cutoff at the edges, and when the viewing direction is offset from the jet, the T factor should be replaced 
by a more general Doppler factor (e.g. Woods &: Loeb 1999; Nakamura 1999; Salmonson 2001). However, if the GRB jets are structured with 
smooth variation of energy per solid angle in different directions (e.g. Rossi, Lazzati & Rees 2002; Zhang & Meszaros 2002; see also Woosley, 
Zhang & Heger 2002, within the context of the collapsar models), this expression has already taken into account the viewing effect, since what 
matter are simply the relevant parameters (e.g. T, L, etc.) along the line of sight (Zhang & Meszaros 2002). 



-1/2 

' (16) 

This is radius-dependent (through different B' at different 

/"2") 1/2 

radii), but in such a form (e.g. 7e oc B 1 ) that the 
comoving typical synchrotron frequency only depends on 
fundamental parameters. In the real problems, which of 
the two typical 7 e 's (je or 7,^) is relevant depends on 
the concrete conditions involved, and on the particle spec- 
trum. Below we will discuss both the low a and high a 
cases for the internal models. 
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2.2.1. Low a: internal shock dissipation 

For the internal shock model, electrons are accelerated 
at the shocks via the Fermi mechanism to achieve a power- 
law distribution with the index —p. The minimum Lorentz 
factor of the injected electrons (i.e. 7! in the above no- 
tation) is 7 e = f(p)(m p /m e )(e e /£ e )9 p ~ SlOe^, where e e 
is the electron equipartition factor which is believed to be 
close to unity (e.g. 0.3) in the internal shocks, £ e is the 
injection fraction of the electrons which could in principle 
be less than unity (Bykov & Meszaros 1996; Daigne & 
Mochkovitch 1998), f(p) = (p- 2)/(p- 1), and the typical 
value 1/6 as p — 2.2 has been adopted to get the numerical 
value. The parameter e x2 = [/(p)/(l/6)](e e /£ e )0p could be 
of order unity, if £ e is below unity for more than an order of 
magnitude. One important feature of the internal shock 
model is that 8 P , and therefore 7 e only depends on the 
relative Lorentz factor between the colliding shells, which 
should not strongly depend on the actual bulk Lorentz 
factor of the shell, T, or the radius of the emission region, 

r oc L 2 . Notice that in this paradigm, ji 2 " 1 (eq.[^6|) defines 
the maximum acceleration energy of the electrons, and is 

much larger than 7,^ (~ 310), and je is relevant in the 
problem. Re- writing (12), we get the peak energy of the 
internal shock synchrotron model, 

E p (low a) ~ 200 keV e^C^Ts'X 1 + Z Y\ (17) 



where e X 3 = e x ie x2 smty . It is worth emphasizing that 
the r _1 dependence in (17) infers that E p oc T~ 2 (see also 



Ramirez-Ruiz & Lloyd- Ronning 2002). 

In some cases, e.g., long t v (large r) or large £ e (so 
that e x2 is much below unity) , the synchrotron peak is be- 
low the BATSE window, and the synchrotron self-inverse 
Compton component may be responsible for the BATSE 
emission (Panaitescu & Meszaros 2000). The peak energy 
in this model is E p c = j 2 E p yn , or 

E p JC {low a) ~ 200 keV e^L^ r^){l + z)" 1 , (18) 

where e X 4 = e x \ (e x2 /0. 06) 4 sin In both eqs.(fL7|) and 
(|T8]), all the parameters related to the shock physics are 
absorbed into two single parameters, e^, e x 4, and the de- 
pendence, E p oc 1^/2 r Xi (1 + z) -1 is valid for both the the 
synchrotron and IC components. Several caveats ought to 
be made for the IC-dominated spectrum. First, es/e e -C 1 
is generally required to achieve a prominent IC component 
(Zhang & Meszaros 2001 and references therein). In the 
internal shock regions, the effective cb can not be too low 
below unity given a not too small wind component B' w (e.g. 
a > 0.1). Second, an IC-dominant fireball is inefficient 
in converting energy to radiation since multi-IC compo- 
nents at even higher energies consume much more energy 
so that one requires a very large energy reservoir to begin 
with (Derishev et al. 2001). Third, since e x z oc e x2 while 
e X 4 oc e x2 , the IC mechanism tends to produce a wider E p 
distribution than the synchrotron mechanism (see §3 and 
Fig.l for more discussions). 

2.2.2. High a: magnetic dissipation 

If a wind has a high a > a c \ ~ (0.1 — 1), strong in- 
ternal shocks may not occur. Alternatively, intense en- 
ergy dissipation may occur through magnetic reconnec- 
tion and/or plasma instability. For a c \ < a < a c2 ~ 190, 



MHD does not break globally, but local dissipation is pos- 
sible through reconncction. If reconnection occurs below 
the photosphere, the dissipated energy is used to accel- 
erate the fireball bulk. Only when a is above some value 
could the energy dissipation occur beyond the photosphere 
(Drenkhahn 2002). For a > a c2 , the MHD approximation 
breaks down globally before the deceleration radius, and 
a more violent, global energy dissipation is possible. The 
particle acceleration mechanism in both cases is hard to 
delineate. For the local dissipation case, the in-principle- 
achievable electron energy may be defined by the DC com- 
ponent of the electric field, i.e., j e ,M,DC m eC 2 = eE'l', 
where I' ~ Tc/fl ~ 10 9 r2^7 is typical comoving scale of the 
stripped wind which may be regarded as the characteristic 
reconnection scale, E' — (x(va/c)B' is the comoving DC 
electric field, where va ~ c is the Alfven speed, and a is a 
poorly constrained parameter (Syrovatskii, 1981; Craig & 
Litvinenko 2002). This gives 7 e ,M,r>c ~ 6 x 10 11 a for typ- 
ical values, which is extremely high unless a is very small. 
In reality, such a linear acceleration would be disrupted by 
any perturbation perpendicular to the DC direction, and 
the maximum achievable j e may be still determined by 
the unavoidable perpendicular magnetic fields as well as 
the radiation self-reaction. Estimating the limiting 7,= is 
difficult due to many uncertainties involved. In any case, 
if the electric and magnetic fields in the dissipation region 
could be regarded as quasi-random, the situation would 
be quite similar to the case of the global magnetic field 
dissipation, which we discuss below. 

For a > o~ c2 , MHD breaks down at r MHD < rd ec , and 
magnetic dissipation occurs globally. One possible conse- 
quence is that the MHD wind is converted to a large ampli- 
tude electromagnetic wave (L AEMW) . The electrons may 
then surf and get accelerated in such a LAEMW (Gunn 
& Ostriker 1971; Michel 1984; Usov 1994; Chen, Tajima 
& Takahashi 2002). Alternatively and more plausibly, the 
LAEMW will soon evolves non-linearly and induces an 
electromagnetic turbulence through an overturn instabil- 
ity, so that the electromagnetic fields in the dissipation 
site are randomized (Lyutikov & Blackman 2001). Parti- 
cles are stochastically accelerated and emit in the random 
fields. To estimate the maximum electron energy gained in 
the random field, we can use the two constraints that lead 

to the estimates 7e and 7^ , respectively. The first con- 
straint is energy conservation. By assuming that a fraction 
e e i of the co- moving local magnetic field energy density is 
eventually converted to the co-moving kinetic energy of 
the electrons, one can always write 



j e m e c 2 ^ e (n' b + n'±) = e e /B' /8n, 



(19) 



where n' b and rij_ are the comoving baryon (and hence, 
baryon associated electron) and pair number densities, re- 
spectively. As discussed in §2.1, for the typical a of inter- 
est, n' b S> n' ± . Without generating additional pairs (which 
seem to be likely as discussed below), the emitting elec- 
trons are only those associated with the baryons, and this 
finally gives 7! ~ 2.8 x 10 5 e xl e e >a 2 ^~ 1 , which has no 
apparent r- or T-dependences since v! oc B oc {L'/r 2 ). 
Comparing to the typical ~ 310 in the internal shock mod- 
els, this 7e is too large to reproduce the observed E p in the 
synchrotron model, unless e e i ~ 10~ 3 <C 1. 
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The second constraint due to radiation self-reaction 
would be more relevant in this case, since the electrons 
are directly accelerated in the random electric magnetic 
fields. A similar estimate as ( |l6| ) can be made. In 
a random electromagnetic field, the acceleration due to 
the electric field parallel to the moving direction (En) is 
much smaller than that due to the electric field perpen- 
dicular to the moving direction (E±) so that the whole 
field may be regarded as a random magnetic field with 
effective strength & ~ \[2B' (Landau & Lifshitz 1975; 
Lyutikov & Blackman 2001). Substituting B' by B' in 
dig), the maximum achievable electron Lorentz factor is 

7e 2 ' ) ~ 3.9 x 10 4 £?' 6 X ^ 2 . This gives the typical comoving 
synchrotron frequency E' syn ~ 30 MeV. Obviously, this 
is too high to meet the observations. 

There are two possibilities. First, it could be that the 
turbulence scale is much less than the radiative scale ct' c /T, 
so that the electrons never reach the maximum accelera- 
tion, but are limited by the turbulence scale itself. This ef- 
fectively lowers the achievable 7 e by a constant factor with- 
out modifying the dependence laws (Lyutikov & Blackman 
2001). One may assume 

le = tc~l¥ ] = 3.9 x 10 2 e c ,_ 2 B' 6 ~ 1/2 , (20) 

so that 

E P (high a) ~ 300 keVe^Cl + (21) 

where e x § — e Ci _2sin , I'. The drawback of this possibility 
is that the process is still inefficient, since a low e c ~ 10 -2 
again implies a very low e e < ~ 1% given n'± -C n' b . This 
is in contrast to the motivation of the magnetic dissipa- 
tion model (e.g. Drenkhahn & Spruit 2002) which aims 
to overcome the low radiation efficiency encountered in 
the internal shock models. The main reason for the low 
efficiency in this case is that for a higher a wind, there 
are less baryon-associated electrons emitting given a same 
total wind luminosity. 

The second (more plausible) possibility is that the pri- 

(2) 

mary electrons do get accelerated close to and emit the 
~ 30 MeV photons in the co-moving frame. These pho- 
tons will interact with the low energy photons (comoving 
energy ~ 10 keV) to produce pairs if the radiation den- 
sity is high enough. The new pairs are accelerated via 
the same mechanism, and a pair cascade develops, un- 
til eventually the pair density n± = £,± n b ^ n 'bi which 

(2) 

limits the achievable 7 e to be below 7e through the en- 
ergy budget constraint Since there are more leptons 
emitting (most being pairs), a higher radiation efficiency 
is achievable. Though the cascade process is hard to de- 
scribe analytically, one may estimate the minimum multi- 
plicity by demanding the typical comoving emission energy 
to be below the pair threshold, i.e., j^(heB' /m e c) <0.5 

MeV. This gives 7e = (0.5/30) 1 / 2 7? ) - 5.0 x 10 3 B' Q 1/2 . 
Assuming a high radiation efficiency, i.e., e e i ~ (1/3) in 
cq. ([l9|) , the required minimum pair multiplicity is £± = 

19cr2e e ',i/3^52 4r i3 1 ^ 2 r2 1 ^ 2 - This still results in too high 
an E p . In reality, more pairs would be further produced 
as long as there is a substantial hard energy spectrum ex- 
tended above the typical synchrotron energy, which would 
degrade the expected E p further. Notice that as long as 



optically thin to the pairs, the dependence 7 e cx B' is 
still valid, mainly because both the comoving pair thresh- 
old itself (~ 0.511 MeV) and the comoving characteristic 
synchrotron frequency (~ 30 MeV for the first generation) 
are constant. The expression ( pl| ) is still valid, with e x $ 
being defined differently. The attractive feature of such a 
model is that the final E p scatter is only due to the dis- 
persion of the bulk Lorentz factor T and a parameter e x §. 
This is favorable for the observed narrow E p distribution 
for the bright BATSE bursts (see §3 and Fig.l for more 
discussions). 

In the magnetic dissipation model, there could be also a 
self-IC component. However, since the synchrotron com- 
ponent already has too high an energy for the GRB emis- 
sion, invoking IC as the GRB mechanism has to introduce 
some ad hoc assumptions. Also all the criticisms to the IC 
internal shock models apply here as well. We therefore do 
not discuss such a possibility further. 

2.2.3. Optically-thick pair photosphere 

For both the shock dissipation and the magnetic dissipa- 
tion scenarios discussed above, there could be the possibil- 
ity that secondary pairs are so abundantly generated that 
they form an optical thick screen for the emergence of the 
photons (e.g. Guetta, Spada & Waxman 2001; Kobayashi, 
Ryde & MacFadyen 2002; Meszaros et al. 2002). There- 
fore the observed GRB emission, at least for some cases, 
has undergone Compton multi-scattering. The emergent 
GRB spectrum is then expected to be regulated by the 
existence of the pairs. The existence of the pairs may 
smear out or even destroy the "clean" correlations dis- 
cussed above. Due to the non-linear nature of the prob- 
lem, simple analytic E p dependences (i.e., the substitution 
of eqs. £3 and pi) are hard to provide without numerical 
simulations. In any case, Meszaros et al. (2002) have 
shown that pairs tend to be self-regulated at a moderate 
optical depth of t' ± ~ a few, and the co-moving typical 

frequency E' cx Tj/t' 2 ^. Since the comoving effective pair 
temperature is insensitive to the shock radius as long as 
copious pairs are produced, such a case is analogous to the 
high-er magnetic dissipation case, but due to quite differ- 
ent reasons. The dependence E p cx T(l + z) _1 is more 
or less retained, with the possible weak dependence on r 
through T± . More detailed numerical work is underway to 
test such a simple treatment (Ramirez- Ruiz et al. 2002, 
in preparation). 

The condition for forming such a pair photosphere is 
subject to further investigation. Pilla & Loeb (1998) nu- 
merically simulated the pair production process in the in- 
ternal shocks with certain model parameters, and found 
that the pair processes mainly distort the high energy part 
of the spectrum, with the low-energy synchrotron peak 
almost unaltered. This hints that at least for some pa- 
rameter regimes (which may be large), the optically-thin 
synchrotron model (both low-cr and high-a cases as dis- 
cussed in §2.2.1, and §2.2.2) applies. Another caveat is 
that a pair-regulated emission component tends to have 
smoother lightcurves, since short time scale variability will 
be smeared out through multi-Compton-scattering (see, 
e.g. Lazzati 2002). This hints that at least for those skipy 
lightcurves with short time variability, such an emission 
component is not important. 
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2.3. External models 

The physical conditions at the deceleration radius, in 
certain circumstances, also allows emission with typical 
energy in the BATSE band, both for the low a and the 
high a cases. The main difference for both cases is the ori- 
gin of the magnetic fields in the energy dissipation regions. 

For a low-cr flow (e.g. a < O.f), the picture is the 
familiar external shock model (Meszaros & Rees 1993). 
The comoving magnetic field is generated in-situ through 
turbulent motions to some fraction of the equipartition 
value, B' = (32Tre B n ext m p c 2 T 2 ) 1 / 2 ~ 40 G e 1 J 2 n 1 J 2 t T 2 . 
The typical electron energy is j e = [(p — 2)/(p — 
l)](m p /m e )(e e /£ e )r = 3.1 x 10 4 (e e /^)[/(p)/(l/6)]r 2 , 
where T is the Lorentz factor of the forward shocked re- 
gion. The E p prediction in this model is 



B' w ~ Q.le Wt -\B' , and that the typical electron energy is 
7 e ~ 200e 7 r ~ 2 x 10 4 e 7 r 2 . This gives the right E p energy 



E e p (low a) = 75 keV e x6 nJ x 2 T A 2 {l + z) 



(22) 



where e x6 = e B 7 (e e /£ e ) 2 [/(p)/(l/6)] 2 sin It is interest- 
ing to note that T appears to the fourth-power, which can 
magnify any small dispersion in T and tend to broaden 
the E p distribution (see §3 and Fig.l for more discus- 
sions). Another emission component in the low-cr external 
shock scenario is emission from the reverse shock propa- 
gating into the ejecta. Compared with forward shocked 
region, the reverse shocked region has the same comoving 
magnetic field strength and the bulk Lorentz factor, but 
the typical random Lorentz factor is 7^ = [{p — 2)/(p — 
l)](m p /m e )(e e /^)r r = 3.1 x 10 2 (e e /£ e )[/(p)/(l/6)]r'', 
where r r ~ Fo/2r is the Lorentz factor of the reverse 
shock in the rest frame of the undecelerated ejecta with 
the bulk Lorentz of r (see e.g. Meszaros & Rees 1997a; 
Sari & Piran 1999). The typical synchrotron frequency 
is therefore degraded with respect to that in the forward 
shock by a factor of T 2 , and lies in the optical band for typ- 
ical parameters. However, the IC emission off these soft 
photon fields by the electrons both in the reverse shock 
and the forward shock regions can result in MeV photons 
which may (partially) account for the GRB prompt emis- 
sion (see Wang, Dai & Lu 2001 for detailed discussions). 
In any case, these IC components have high powers of the 
random Lorentz factors (either cx r 2 (r r ) 4 or cx r 4 (r r ) 2 ), 
and will have a similar or even broader En distribution as 
the standard external shock model (eq. |22| ) . Hereafter we 
do not explicitly discuss these possibilities. 

For the high-cr case, the external shock variant is given 
by the interaction between a relativistic plasma stream 
and a magnetic barrier. As simulated by Smolsky & Usov 
(2000), the outcome is dependent on the energy density ra- 
tio between the plasma and the field, a = 8Trn ex tm p c 2 (T — 
1)/B' 2 . In such a case, the comoving magnetic field B' is 
dominated by the primordial wind component at the de- 
celeration radius, i.e., B' ~ 180 G L^^E^^T^^n],^. 
The energy of electrons accelerated in the electric field 
generated at the front of the magnetic barrier is typi- 
cally 7 e ~ (m p /m e )T — 1.8 x 10 5 r 2 . The typical syn- 
chrotron energy is too high, and Smolsky & Usov (2000) 
invoke the synchro-Compton radiation in a LAEMW as 
the GRB radiation mechanism. The latter is essentially 
a synchrotron mechanism except that the comoving B' is 
now replaced as the amplitude of the LAEMW, which is 



E p (high a) 



880 keV e x7 r 8 2 /3 L^E 5 ^*n^ t (l + z) 



1/3/ 



-V2 P -V3„ 
^52 -^52 

(23) 

where e x i = e w .^\e~,a^ 2 {\ + cr) 1 / 6 . The main difference 
between (^) and ( [22] ) is the origin of B'. The wind B' is 
determined by the properties of the central engine, while 
the shock generated B' depends on T. Thus the power of 
r in (23) is 4/3 lower than that in (22), leading to less 
dispersion in E p for the same T scatter (see Fig.l). 

In both cases, pair production is likely not to be crucial 
in regulating the E p , due to a much smaller compactness 
parameter involved. 

2.4. Innermost models 

Regardless of the a value, a prompt signal is emitted as 
soon as the fireball becomes Thomson thin. The typical 
energy of such a baryonic photosphere emission is defined 
by the photosphere temperature, which depends on the 
dimcnsionless entropy of the fireball. Our treatment be- 
low follows Meszaros et al. (2002), with the modification 
that we include a Poynting flux component. Since this 
cold component is not in the form of photons initially, it 
is left out in calculating the temperature of the baryonic 
photosphere. 

The compactness parameter of the fireball can be 
expressed as a dimensionless parameter 6 , l Vj0 = 
(LaT/^rn p c 3 r ) — 1.2 x 10 12 L5 2 r^~7. For the pur- 
pose of the following discussion, L should be substi- 
tuted by the "hot" component, L/(l + a). The photo- 
sphere radius is defined by the Thomson-thin condition, 
i.e., UTfi'V = 1, where I' is the comoving length of the 
continuous wind or the discrete shell. There are two 
critical values of the dimensionless entropy in the prob- 
lem, both of which are related to £ Pt0 . The first value 

T} cl = el{o = 2.5 x 10 2 [L 52 (1 + er)- 1 ^] 1 / 5 is the criti- 
cal entropy below which the opacity is defined by a con- 
tinuous wind rather than a discrete shell. The second 
value r? c2 = il(o = 1.0 x 10 4 [L52(1 + v)' 1 ^}] 1 / 3 is the 
critical entropy above which the photosphere occurs in 
the acceleration regime. The observer-frame photosphere 

temperature satisfies p h/Oo = (t7 8 / 3 77~ l 10/ ' 3 , 77/^2, 1) for 
(77 < 77 c i,77 c i < i] < rj c2 ,ri > 77 c2 ), respectively (eq.[6] in 
Meszaros et al. 2002, see detailed derivation given there). 
If the GRB emission is due to Comptonization of the pho- 
tospheric emission by the Alfven turbulence (Thompson 
1994), the resultant E p should resemble the thermal peak 
predicted by Q p h- Noticing (Q), we get 

E p (ph)(l + z) = 2.8Q ph m e c 2 = 
11 keV L- 2 5/12 (l + cr) 5 / 12 ^l„ 3 r« /3 , V < Vcl] 
370 keV L- 2 1/12 (l + a^H-^U^, Vd < V < vd^) 
2.7 MeV Ll{\l + 0)^/%^%, V > Vc2 . 

This is an extension of the discussion in Thompson 
(1994), who only discussed the highest rj case. The to- 
tal photosphere luminosity scales with the temperature, 



3 Conventionally, m e rather than m p is used to define I, e.g. Guibert, Fabian & Rees (1983). 
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i.e. Lph/Lo = &ph/Qo- Its relative importance with re- 
spect to the internal component depends on the regimes 
where r\ lies. For the low-cr case, the internal shock com- 
ponent dominates when r\ is small, e.g. rj < r] c x, but 
the photosphere component starts to be prominent when 

7? >^/4^1 3(l + cr )-V3. 

3. TESTING MODELS AGAINST DATA 

We consider the following £" p -related observational 
properties: 

(i) For the 5500 spectra of 156 bright BATSE bursts an- 
alyzed in Preece et al. (2000), the typical E p is 300 keV, 
and the distribution is a narrow lognormal with full-width 
at half-maximum of less than a decade. Since each burst 
has many spectra in this dataset, the result implies that 
the E p dispersion for a particular burst is also very narrow. 

(ii) A substantial population of X-ray flashes has been 
recently studied, with typical E p < 40 keV (Heise et al. 
2001; Kippen et al. 2001). If these events have the same 
origin as normal GRBs, they extend the E p distribution 
to a wider range, but the small sample collected so far 
indicates that the log-normal shape seem to be violated 
(Kippen et al. 2002). 

(iii) Statistically, a positive correlation between E p and 
the isotropic luminosity L, E p oc L s with 5 ~ 0.5, has been 
noted recently (Amati et al. for the BeppoSAX bursts, 
and Lloyd-Ronning & Ramirez-Ruiz 2002 for the BATSE 
bursts). 

These properties, along with other information such as 
temporal variability, spectral lags, etc. are relevant quan- 
tities to be accounted for in a theoretical model. Below 
we will test the model E p predictions with the above ob- 
servational facts. The model predictions are collected in 
Table §. 

3.1. Narrowness of E p distribution 

To compare the narrowness of the E p distributions 
among different models, we performed a simple Monte 
Carlo simulation. Since we are concerned about the nar- 
rowness of the distribution rather than the absolute value 
of E p in different models (and the uncertain free parame- 
ters make it possible for all the models to adjust the abso- 
lute E p values to match the 300 keV typical value), we have 
adopted arbitrary units for the parameters so that log E p 
peaks around 0. We randomly generated 40000 bursts with 
distributions of the unknown parameters as follows: (1) 
burst luminosities span three orders of magnitude with a 
N(L)dL oc L~ 2 dL distribution (such a luminosity func- 
tion is the natural consequence of the universal jet model, 
Zhang & Meszaros 2002, Rossi et al. 2002; and has re- 
ceived preliminary support from the data, Schaefer, Deng 
& Band 2001; Schmidt 2001); (2) burst Lorentz factors 
span 1.5 orders of magnitude and are correlated with the 
luminosity as T oc L 1 / 2 with random deviations (a positive 
r — L correlation is expected in both the conventional and 
the universal jet models, e.g. Kobayashi et al. 2002; Rossi 
et al. 2002; Zhang & Meszaros 2002, and here we adopt 
the (1/2) power as an example); (3) burst redshifts are ran- 
domly distributed in the (0 — 15) range (a more detailed 
distribution function following the star formation history 
of the universe essentially does not influence the width of 
the final E p distribution, introducing only a minor dis- 



tortion of the distribution shape); (4) logi„, log(T4/rj_ 2 ), 
l°g e a;5, log£a;6, log £x7j are randomly distributed between 
(0,1); (5) \ogn ext is randomly distributed between (0,2), 
and the total duration logT is randomly distributed be- 
tween (0,3); (6) foge^i and log e X 2 follow a standard Gaus- 
sian distribution with random deviations. (This is because 
the product of at least three random variables has a log- 
normal distribution, Ioka & Nakamura 2002). For com- 
parison, we also let e X 3 = e x ie X 2 be represented by a unit 
lognormal itself for another realization. The E p predic- 
tions within the different models are calculated according 
to Table |[ The fluxes of the bursts are calculated ac- 
cording to the standard (Cl m )&A) = (0.3,0.7) cosmology. 
Given a certain flux threshold, the histograms of the E p 
distributions for different models are plotted for compari- 
son. 

Figure [l] shows the E p distributions of various models. 
To take into account the flux threshold effect, only the 
bursts whose fluxes are among the top-three flux decades 
are selected (2962 bursts from a total of 40000 bursts sim- 
ulated). For ease of comparison, we have plotted all the 
histograms with the same scale (which spans 10 orders 
of magnitude). The models in the different subplots are: 
(a) the internal shock synchrotron model, with both e x i 
and e X 2 distributed as a lognormal; (b) the same inter- 
nal shock synchrotron model, but with e x s itself assumed 
to be a standard lognormal. (This would imply that the 
dispersions among 9 p , e e , es, etc. are small due to some 
unknown conspiracy); (c) the internal shock synchrotron- 
self-IC model corresponding to the synchrotron model (a); 
(d) the internal magnetic dissipation model, which also ap- 
plies to the pair-dominated internal models; (e) external 
shock model; (f) external magnetic model; (g) baryonic 
photosphere model in the coasting wind regime; (i) bary- 
onic photosphere model in the coasting shell regime; (j) 
baryonic photosphere model in the acceleration regime. 

From Fig.[j] the following conclusions can be drawn. 1. 
Given the assumed parameter distributions, neither the 
internal nor the external models can reproduce the nar- 
row E p distribution found in the sample of bright BATSE 
bursts (Preece et al. 2000). In our simulations, it is 
found that the narrowness of the E p distribution depends 
on the number of the independent parameters involved 
in the model, the narrowness of the assumed distribution 
of each parameter, as well as the power of these param- 
eters. The more independent parameters, or the higher 
the power to which a certain parameter is involved in a 
model, the broader is the resulting E p distribution. 2. 
The IC models (Fig.lc) and the external models (Fig.le 
and If) are even less favored, due to the much wider dis- 
tributions caused by high powers of some parameters in 
their E p expressions (e.g. oc e x2 for the internal shock IC 
model, and oc T 4 in the external shock model). 3. The 
usual internal shock synchrotron model (Fig. la) also in- 
volves many parameters, and therefore results in a broad 
E p distribution among bursts. However, assuming that 
some of the parameters are not independent, one can re- 
duce the number of free parameters to obtain a narrower 
distribution (Fig. lb). 4. Both the internal magnetic model 
and the pair-dominated model invoke the least indepen- 
dent parameters, and have narrow E p distributions. How- 
ever the simulated distributions are still about one order 
of magnitude broader than the E p distribution observed 
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in the bright BATSE bursts. 5. Two sub-classes of the 
innermost model, (h) and (i), have very narrow E p dis- 
tribution (due to much smaller powers of the parameters 
involved, eq.[Q). These are comparable to that of the 
bright BATSE bursts. However, it is likely that these 
emission components contribute only partially, and under 
certain circumstances, to the GRB prompt emission 

Several caveats about our simulations are pertinent. 
First, our primary goal here is not trying to match the 
Preece et al results closely. To make such close compar- 
isons, other effects such as the detector energy band and 
the instrument response function ought to be taken into 
account, and a much larger parameter space would need 
to be explored for every model. It is also possible that 
the intrinsic E p distribution is broader than what is found 
by Preece et al. (2000), as indicated by the growing pop- 
ulation of the XRF sources. Instead, we have adopted a 
uniform set of parameters for different models, which al- 
lows an unbiased evaluation among the models. In fact, it 
may be possible to reproduce Preece et al's results by ad- 
justing the distribution function and the dispersion of the 
input parameters within a certain model. For example, 
Bottcher & Dermer (2000) have shown that Preece et al's 
narrow E p distribution is reproducible even in the E p oc T 4 
external shock model, by taking into account the instru- 
mental effects and the flux threshold effect, although they 
had to fix some input parameters (e.g. n ext ). Our Figfil 
nonetheless clearly reveals the "relative" abilities of differ- 
ent models to generate a narrow E p distribution. Second, 
in our simulations most of the parameters are regarded 
as independent of each other. However, there may exist 
some intrinsic correlations among some parameters, and 
the conspiracy among these parameters can reduce the E p 
scatter to achieve narrower distributions (we have shown 
this effect in Fig. la and lb). More detailed work is needed 
to explore these intrinsic correlations. Finally, it is possi- 
ble that the flux threshold influences the scatter of the dis- 
tribution. In Figj2[ we have explored such flux threshold 
effects for a particular model (the internal magnetic model, 
same conclusion also applies to other models). The three 
subplots indicate different flux thresholds: (a) E p distri- 
bution of all the 40000 bursts; (b) that of the 2962 bursts 
whose fluxes are among the top three orders of magnitude; 
(c) that of the 848 bursts whose fluxes are among the top 
two orders of magnitude. We see that the E p distribution 
indeed gets narrower when raising the flux threshold, but 
the effect is only mild. This might be because the intrinsic 
E p scatter due to many unknown parameters has a similar 
effect in a very wide flux decades. 

3.2. Nature of X-ray flashes 

Following the reports of the identification of the XRFs as 
a type of cosmological explosions resembling GRBs (Hcisc 
et al. 2001, 2002; Kippen et al. 2001; 2002), several specu- 
lations on the nature of these objects have been proposed. 
Whether they smoothly join the GRB distribution or con- 
stitute a related but separate class, the E p distributions 
associated with these objects would depend on the possible 
models, some of which we discuss below. 

1. Dirty fireballs? Based on the external shock model, 
Dermer et al. (1999) discussed a type of then-undiscovered 
objects with lower bulk Lorentz factors than conventional 
GRBs, and suggested that such dirty fireballs would typi- 



cally have lower E p (due to the E p oc T 4 dependence) . It is 
then natural to attribute the XRFs to such dirty fireballs 
(Heise et al. 2001). However, we note that a dirty fireball 
does not necessarily produce a low E p burst. In the opti- 
cally thin internal shock model, equation ( |l5| ) shows that 
the apparent T dependence cancels out, and equation ( p~T[ ) 
indicates that a dirty fireball has a closer-in internal shock 
radius, hence a higher E p due to a higher magnetic field 
(Fig. 3a). In the high a case or the pair-dominated case, 
and in the external models, a positive dependence of E p 
on r is retained, and a dirtier fireball tends to produce a 
lower E p (Fig.|). 

2. High redshifl bursts? An interesting possibility is that 
XRFs are GRBs at much higher redshifts (e.g. z > 6), 
which may thus be related to the death of the first stars 
in the universe (Heise et al. 2001). However, some of the 
other collateral evidence for high-z location, e.g., time dila- 
tion for both the total durations and the individual pulses, 
is lacking, which casts some doubts on such an interpre- 
tation (Lloyd- Ronning 2002). This problem is made less 
severe by the fact that also in normal GRBs, the presence 
of such time dilation is less noticeable than the intrinsic 
dispersion, and becomes noticeable only when analyzing 
very large samples (e.g. Norris et al. 2000). 

3. Off-beam GRBs? There is now strong evidence that 
at least a fair fraction of GRBs are collimated. This raises 
the possibility of interpreting the XRF phenomenon as be- 
ing due to viewing angle effects. In a version of such mod- 
els where a sharp jet edge is assumed (Yamazaki, Ioka & 
Nakamura 2002), XRFs are those off-beam GRBs in which 
the observer misses the bright jet cone. Such a model re- 
quires the XRFs to be very nearby z<0.2. The redshift 
measurement of one X-ray rich GRB 011211 (a close rel- 
ative of XRFs) at z = 2.14 seems to suggest that at least 
some XRFs may not be interpreted in this scenario. In 
the universal jet beam model (Rossi et al. 2002; Zhang 
& Meszaros 2002), the off-beam population is greatly re- 
duced, but bursts observed at large viewing angles tend 
to be "dirty". The possible connection of such a config- 
uration to XRFs (within the framework of the collapsar 
models) has been suggested by Woosley et al. (2002 and 
references therein). 

4. Photosphere- dominated fireballs? Meszaros et al. 
(2002) suggest that XRFs may be accounted for within the 
standard fireball internal shock model with moderate red- 
shifts, due to a dominant contribution of either the bary- 
onic or the shock pair photosphere. A similar proposal in 
the high-(7 case has also been suggested (Drenkhahn 2002). 

The first two speculations invoke the dispersion of one 
particular parameter, while the total E p dispersion is 
the combined dispersions of many independent parame- 
ters. Whether the influence of one particular parameter is 
prominent depends on how many independent parameters 
are involved, and what is the power of that particular pa- 
rameter in the E p expression. To show this effect, we plot 
in FigJ| the predicted E p as a function of T in four models: 
(a) internal shock model; (b) internal magnetic model; (c) 
external shock model; and (d) external magnetic model. 
We find that the E p broadening due to the T scatter is 
more evident for the external models in which high powers 
of r are invoked. Such an effect also depends on the band- 
width of the detector. Only when the bandwidth is larger 
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than the intrinsic scatter at a same T could the broadening 
due to r scatter be noticeable. For the same reason, the 
redshift scatter only contributes weakly to the E p scat- 
ter, so that the high redshift interpretation of XRFs is not 
quite plausible. The photosphere interpretation in fact in- 
vokes possibly two (or even more) emission components 
with different E p distributions. A natural expectation of 
this model is that the E p distribution of the GRB/XRF 
combined population may have a non-lognormal or even 
bimodal distribution, which seems not incompatible with 
the current small sample (Kippen et al. 2002). To iden- 
tify the nature of the XRFs, clearly more data are needed, 
including both redshift measurements and spectral fits. 

3.3. E p — L correlation 

The positive correlation E p cx L s , (Amati et al. 2002; 
Lloyd-Ronning & Ramirez-Ruiz 2002) also poses interest- 
ing constraints on the models. It is seen in Table |^ that the 
different models give different predictions for the E p — L 
dependence. The complications arise from the E p depen- 
dence on the unknown T expected in these various models, 
since there is no simple way to relate T to the observables. 
In the above Monte Carlo simulations, we have adopted a 
positive r cx L k (k — 1/2 specifically) correlation, which 
is consistent with the expectation of the theoretical mod- 
els (Kobayashi et al. 2002; Rossi et al. 2002; Zhang & 
Meszaros 2002; Salmonson & Galama 2002). If one takes 
k as a free parameter, some constraints on k may be im- 
posed by the observed S. For example, the internal shock 
synchrotron model could not give a positive E p — L corre- 
lation if k > 1/4. Ramirez- Ruiz & Lloyd-Ronning (2002) 
also noticed this, and argued that their observed rela- 
tion (Lloyd-Ronning & Ramirez-Ruiz 2002) can be only 
accommodated within the IC-dominated model with an 
additional assumption of a 7 e cx T (the T~ 2 dependence 
thereby canceling out). However, we note that the IC- 
dominant picture is less favored for other reasons as dis- 
cussed above. Similar analyses could be made for the other 
models. For example, both external models seem to pre- 
dict a too steep E p — L correlation, unless the T—L relation 
is rather mild. The internal magnetic model or the pair 
dominated model are however compatible with the data 
for reasonable parameters. 

4. CONCLUSIONS 

The nature of the GRB prompt 7-ray emission, includ- 
ing the emission site and the energy contents of the fire- 
ball, are still poorly known after more than 30 years of ef- 
forts. We have analyzed the various fireball model variants 
within a unified picture, and have revisited the E p predic- 
tions of different models. These models are tested against 
the current GRB spectral data with a simple Monte Carlo 
simulation, with attention to the narrowness of the distri- 
bution and its dependence on some of the physical param- 
eters. Our aim is to set up a general theoretical framework 
to allow unbiased tests of these models against the known 
data. Based on the analysis in this paper, we can evaluate 
the existing fireball models as follows. 

1. The internal shock model is generally regarded as 
the most attractive candidate for the prompt 7-ray emis- 

7 It is worth mentioning that within the same burst, the internal 
incompatible with Preece et al. (2000) unless some fine tuning is made 



sion of classical GRBs. The highly variable, spiky GRB 
lightcurves are naturally reproduced in such a model, and 
many studies have shown that it is successful in repro- 
ducing many of the GRB properties (e.g. Kobayashi, Pi- 
ran & Sari 1997; Daigne & Mochkovitch 1998; Panaitescu, 
Spada, Meszaros 1999; Spada, Panaitescu, Meszaros 2000; 
Ramirez-Ruiz & Fenimore 2000; Guetta et al. 2001). Our 
findings in this paper indicate two important caveats for 
the internal shock model. First, unless the dispersions in 
the shock parameters (e.g. e e , es, 9 P , p, etc) are very 
small or there exist some intrinsic correlations among the 
parameters, the internal shock model generates an E p dis- 
tribution (Fig. la) which, even in the optimistic case, is at 
least one order of magnitude broader than the straight (i.e. 
not otherwise corrected) distribution given by Preece et al. 
(2000). The calculated distribution may be still compati- 
ble with the data if XRFs are included in the observations, 
which intrinsically broaden the E p distribution . Second, 
the synchrotron internal shock model may not be able to 
reproduce the E p — L positive dependence unless the V 
distribution is un-related or weakly related to L, or some 
further assumption is made (e.g. Ramirez-Ruiz & Lloyd- 
Ronning 2002). 

2. The high-cr internal model could inherit most of the 
merits of the internal shock model, with some additional 
advantages such as a smaller E p dispersion and the right 
E p — L correlation. The caveat in such models is that they 
have so far been less well-studied than the internal shock 
models, including the basic particle acceleration and emis- 
sion processes. It is also unclear how to lower the typical 
E p energy to the sub-MeV band, and the physical context 
of how a high-cr flow is launched in a collapsar is not well 
explored. More investigations in this direction are desir- 
able (e.g. Blandford 2002). 

3. In both the shock and magnetic dissipation inter- 
nal models, both a narrow E p distribution and the right 
Ep — L correlation are attainable if a pair photosphere is 
formed. However, it is unclear how common such a situ- 
ation would be. The sharp lightcurve spikes may be also 
hard to reproduce. 

4. The prompt 7-ray emission predicted in external 
models is expected in most cases, unless the external 
medium is very under-dense or previous internal dissi- 
pation of the bulk kinetic energy has been highly effi- 
cient. Whether the observed GRB prompt emission is 
attributable to this component is in question. Our sim- 
ulations show that these models produce too broad an 
Ep distribution compared with other models, and prob- 
ably a too steep E p — L correlation. Other arguments 
against such scenarios include the need for additional as- 
sumptions (blobs) and the inefficiency involved in inter- 
preting the variability (Sari & Piran 1997, but see Dermer 
& Mitman 1999). An important feature of these models 
is that E p are positively dependent on the ambient ISM 
density n ext , which is in principle a measurable parame- 
ter. This provides a potential test for the model. More 
broadband afterglow fits are needed before a significant 
statistical evaluation can be made. Even if the external 
scenario does not account for the GRB prompt emission, 
studies of such a model are nevertheless meaningful since 
they can explore the important bridge between the prompt 

shock model also tends to produce a wide E p distribution, which is 
(S. Kobayashi, 2002, in preparation). 
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emission and the afterglow. 

5. The emission component coming directly from the 
baryon photosphere is expected to partially contribute to 
the observed GRB prompt emission. The E p distributions 
for different regimes are narrow, and E p — L correlations 
are easily accommodated at least in the regime where the 
emission is the strongest (i.e., 77 > 77 c2 ). The best guess is 
that such a component appears mixed in with other com- 
ponents, and becomes important under certain conditions. 

6. All the IC models tend to generate broader E p distri- 
butions as compared with their synchrotron counterparts. 
They are also less favored due to other reasons. 

Finally, our entire discussion in this paper is within the 
context of a naked central engine. A more realistic sce- 



nario invokes the fireball-progenitor envelope interaction, 
e.g. in collapsar models, which is receiving increased atten- 
tion. These would lead to additional emission components, 
which are beyond the scope of the current paper. 
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Table 1 
Fireball variants 



Condition Fireball properties 

<t<o c \ kinetic energy dominated, strong shocks possible 

Cci < a < cr c2 Poynting dominated, MHD does not break globally, but may break locally, no strong shocks 

c C 2 < cr < (T c3 completely Poynting dominated, MHD breaks globally, no strong shocks, baryon electrons non-negligible 

<y > &c3 completely Poynting dominated, MHD breaks globally, no strong shocks, baryon electrons negligible 



Table 2 
Model predictions for E p 



Models 


Sub-categories 


E p 


Comments 


Internal models 


low-cr: internal shocks 
high-c: magnetic dissipation 
pair photosphere 


oce :c3 L 1 /2 r -2 t -i (1 + z) -i 

(xe x5 T(l + z)- 1 
ozTT'Jt'I^ + z)- 1 


7e^ relevant 
7e 2 ' ) relevant 
Comptonizcd spectrum 


External models 


low-u: external shocks 

high-cr: plasma-barrier interaction 


cxe x6 r 4 n^(l + 0)- 1 
oce l7 r 8 / 3 L 1 / 2 £- 1 /3 Il ^(l + z )- 1 


B generated in-situ 

B carried from the wind 


Innermost models 


77 <7j c i~ 250(1 + ( j)- 1 /5 

r/ci < V < Vc2 

f]>Vc2- 10 4 (1 + ct)" 1 ^ 


ocL-^^rS/^l + z)- 1 

ocL-'/^rfiiz)- 1 

ocl^d'tl + z)- 1 


wind coasting regime 
shell coasting regime 
shell acceleration regime 
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Fig. 1. — E p distributions from Monte Carlo simulations for different models: (a) internal shock synchrotron model; (b) refined internal 
shock synchrotron model with a lognormal distribution of e x z\ (c) internal shock IC model; (d) internal magnetic dissipation model, or 
pair-dominated model; (e) external shock model; (f) external magnetic model; (g) baryonic photosphere model in the wind coasting regime; 
(h) baryonic photosphere model in the shell coasting regime; (i) baryonic photosphere model in the acceleration regime. The input parameters 
are described in the text. The histograms represent 2962 bursts whose fluxes are within the top three orders of magnitude, selected from a 
total 40000 bursts. Normalized E p values are used, with typical log E p ~ 0. 




Fig. 2. — The flux-threshold effect on the E p distribution in the simulations for Model (d) in Fig.l. (a) all 40000 bursts simulated; (b) 
the 2962 bursts whose fluxes are within the top three orders of magnitude; (c) the 848 bursts whose fluxes are within the top two orders of 
magnitude. 
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